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■ For the study of reactions in High Energy Physics (HEP) automatic computation systems 

have been developed and are widely used nowadays. GRACE is one of such systems and 
it has achieved much success in analyzing experimental data. Since we deal with the cross 
section whose value can be given by calculating hundreds of Feynman diagrams, we manage 
the large scale calculation, so that effective symbolic manipulation, the treat of singularity in 
the numerical integration are required. The talk will describe the software design of GRACE 
system and computational techniques in the GRACE. 
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Q_i' §1. Introduction 



High energy experimental physics has produced excellent results thanks to the 
progress of the detectors and the development of the high energy accelerators with 
^ ■ high luminosity. Accordingly, the accurate theoretical computation is required to 

compare experimental results with theoretical predictions. On the other hand, as 
the beam energy becomes higher, there appear the physics processes with many 
final state particles. As a matter of course the number of relevant Feynman dia- 
grams becomes huge. This means that Feynman amplitude calculation is practically 
impossible when one calculates cross sections by hands. 

Since the perturbative calculation in quantum field theory is well-defined, it can 
be realized as an automatic computation system. There appeared several systems, for 
instance, CompHEPEP and FeynArt/FeynCalcl$, in HEP and nowadays the automatic 
computation of Feynman amplitudes becomes common. 

We have developed an automatic computation system named GRACE. & In 
section 2, we describe the feature of GRACE system as a concrete example of such 
system. In section 3, several physical achievements by GRACE are briefly explained. 
In section 4, we discuss the techniques in GRACE from the view point of computing. 



* Presented by F.Yuasa at the ICCP5 conference, Kanazawa, October 1999 
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In the final section, we summarize the problems and the future plans. 

§2. GRACE System 

2.1. Tree-Level 

Using GRACE, the computation proceeds in following steps: 

1. define a physics process with a model and generate Feynman diagrams relevant 
to the process, tr 

2. generate FORTRAN code for the numerical computation referring to helicity 
amplitude library (CHANEL) i> 

3. compute the cross section by Monte Carlo integration package (BASES), EP 

4. simulate the events by event generation package (SPRING). B 

GRACE system consists of several component-programs corresponding to each 
step. In the first step program named grc generates diagrams. It can generate di- 
agrams at any order of perturbation. All the resultant information on the diagram 
topology, particles and vertices is stored in the list-formated file. In the next step pro- 
gram named grcf ort generates FORTRAN code including interfaces to CHANEL 
and BASES. In the third step, two programs integ and spring are produced. Be- 
fore proceeding to the third step we have to care about kinematics for the phase 
space integration using BASES and SPRING. Running integ, the cross section is 
computed numerically. Since the distribution of the integrand is known in course of 
the integration, unweighted events are produced by spring in the last step. 

In GRACE, several physics model files and kinematics routines are necessary 
and standard one's are built-in. User-defined model file or user's kinematics file can 
replace built-in files. GRACE also includes utility programs such as diagram-drawer, 
gracef ig, and so on. 

2.2. One-Loop Level 

The important extension of GRACE is to compute the loop amplitudes auto- 
matically. As for the one-loop order, automatic computation of the processes with 
two final state particles (2 — > 2) is almost completed and some 2—5-3 processes 
are computed by GRACE. However, it is currently restricted to the electro-weak 
interaction. 

The steps for the one-loop computation is almost same as those for the tree- 
level. However, extra steps are included due to new features needed for treating loop 
diagrams. We introduced the symbolic manipulation into GRACE to treat them and 
it processes the following steps: 

(i) take the trace if there is a fermion line, 

(ii) introduce Feynman parameters integration method, 

(iii) shift a loop momentum appropriately, 

(iv) drop odd orders of I, 

(v) replace even orders of I as IH V — > g^l 2 jn and so on, 

(vi) contract vertex indices inside the loop in n-dimension. (n = 4 — 2e) 

In order to regularize the ultra-violet divergence and the infrared divergence, the 
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dimensional regularization and the fictitious photon mass parameter A are introduced 
into GRACE, respectively. 

For the loop process program grc generates both relevant tree diagrams and loop 
diagrams at once, grc also generates diagrams with vertex and propagator counter 
terms. After generating diagrams, following steps are needed as a part of step 2. in 
the tree-level procedure: 

• Generate source program for the symbolic manipulation system, REDUCElZP 
or FORMcP. Here, the product r p loop T tree ^ are written down in a symbolic code 
for each pairs of an one-loop amplitude and a tree amplitude. 

• Invoke REDUCE or FORM and get FORTRAN source code. 

After we get FORTRAN source code, the cross sections can be calculated with the 
loop library and the counter term library to deal with Feynman integral and counter 
terms, respectively. 

2.3. How to check the results 

The important and serious issue in automatic computation is how we can check 
the results. For the tree-level computation GRACE can check the gauge invariance 
of the results by using unitary and covariant gauge. This does not work in the loop 
amplitude since the structure of the numerator becomes complex. For one-loop level, 
we newly implemented the non-linear gauge (NLG)cP for gauge invariance check. We 
take a generalized NLG fixing condition for SM: 



{dp - ieaA^ - ig cos 9 w pZ^W +l1 + ^|(« + 5H + ik X3 )x + 

-,}-Uz» + t z -^—{v+~eH) X X --L(^) 2 
2£z V 2costV / 2t4 

(2T) 

Using NLG gauge invariance, check has been performed for several one-loop 
processes. 

Besides gauge invariance check we can check the self-consistency of the loop- 
implementation to GRACE by changing parameters 1/e or A. Comparison among 
the results derived^ from independent automatic computation systems shows very 
good agreements. 

2.4. Beyond the Standard Model 

GRACE has been extended to compute the physics processes not only in the 
Standard Model but also the Minimal Super symmetric extension of the Standard 
Model(MSSM). Here the treatment of the Majorana particles and the fermion num- 
ber clashing vertices is introduced. c3 The MSSM includes 85 particles and 3,764 
vertices so that the automatic calculation becomes indispensable for the calculation 
of SUSY processes. Recently the complete first-order radiative corrections to the pro- 
cess H + -> tb has been calculated using 1-loop system(GRACE/SUSY/lLOOP).B ) 
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§3. Physics Achievements 

So far GRACE system has been proved to be very powerful in computing many 
complicated processes with many final state particles. In the tree-level 2 — ±3 (over 
50 processes), 2 — > 4 (about 100 processes), 2^5 (e + e~ — > e~ i7 e udj i$ ) , and 
2^6 (e + e~ — > bbudixv^ E-^P ) were computed. Among them it is noteworthy that 76 
e + e~ — > 4-fermion processes which take place in the energy region above W + W~ 
pair threshold are computed for LEP-II experiment. These 76 4-fermion programs 
are combined and made them up to the program named grc4f . t3 The one-loop level 
processes such as e + e" —> Z°H, e + e~ —> tt, 77 — > W + W~, and e + e~ — ► W + fiv^ 
were computed. US 

To provide a practical generator for SUSY processes at LEP-II and a linear 
collider, we developed susy23 which contains 23 processes for e + e~ physics. c3 
GRACE/SUSY can be extended to include R-parityviolating interactions so that 
the process e + e~ — ► vv e dd in LEP-II was calculated.!^ 

Recently GRACE has been applied to a process at HERA experiments(GRAPE). 
It was the lepton pair production, both in the elastic and the deep inelastic 
scattering regions, ep — > epZZ and — > eiZX. Some routines for the proton form 
factors and the structure functions have been implemented. 

§4. Issues Related to Computation 

The first issue is the computation time. As the physics process treated by 
GRACE becomes more complicated, the computation time becomes so longer that it 
restricts applications of the automatic computation. Among the computing steps, the 
multi-dimensional Monte Carlo integration consumes a lot of computing time. The 
requirements for reducing the execution time, we have successfully developed the first 
version of parallelized GRACE by using message passing library, PVM or A/TPJ P°P 
This version of parallelized GRACE reduced the computation time drastically. c3 

We measured the reduction rate with Fujitsu AP3000u by using up to 16 pro- 
cessors for the physics process e + e~ — > bbudu^fj,. The speed up with 16 processors 
was 15.82. 

The second issue is related to a treatment of singularities appearing in the multi- 
dimensional Monte Carlo integration. In the framework of GRACE kinematics rou- 
tine used in the phase space integration is not generated automatically. The aim of 
kinematics routine is to find an appropriate transformation between integral vari- 
ables and kinematics variables to keep away from singular behavior of amplitudes 
in the integral region. However, we don't know yet an algorithm good enough to 
deal with general singularities. To find better algorithm for singular integrals is thus 
another big issue and it belongs to the field of numerical integration on computer. c3 

The third issue is the symbolic manipulation. 
• When we write REDUCE code naively, the short of memory often happens dur- 



* The AP3000 system consists of UltraSPARC -II 300MHz processors connected via AP-Net. 
AP-Net is a two-dimensional torus network and provides 200MB/s bandwidth per port. 
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ing processing the manipulation. In order to get rid of this, many complicated 
and sophisticated techniques have been introduced. However in manipulation 
for 1-loop of 2 — > 4 process, it is almost impossible to perform the symbolic ma- 
nipulation by REDUCE. FORM system can handle the lengthy code because 
FORM needs less memory and computing time than REDUCE. 

• When we perform the symbolic manipulation using REDUCE and FORM, the 
long expressions (more than 10,000 lines) are obtained. Then very long comput- 
ing time is required in numerical integration. The effective optimization should 
be considered. The symbolical optimization can be possible through pattern 
matching because Feynman diagrams contain many common parts. 

The last issue is the high precision calculation in GRACE. It is sometime required 
from the following reasons: 

• The all vertices of SM and MSSM had been checked by gauge invar iance in 
the quadruple precision. A possible bug might happen in the kinematics code. 
It can be detected by the gauge invariance check. However, if it affects minor 
contribution, the double precision is not enough to manifest the bug. 

• Due to several kinematical reasons the accuracy of numerical values may be 
lost. For instance, the two photon process has drastic numerical cancellation 
between Feynman diagrams. It may give a wrong result if one uses the double 
precision (10 -15 ). 

• In the cross section formula, there can be a term (W/m e ) 2 which is ~ 10 +14 
in TeV collisions where m e is the electron mass and W is CM. energy. This 
indicates that double precision is no more enough in HEP computation. In 
loop calculation one will loose some accuracy in the evaluation of loop integrals 
in addition to kinematics. In some case the calculation with full quadruple 
precision will be required. 

§5. Summary 

We have developed GRACE, an automatic computation system in HEP, and it 
has made much success for the study of physics process in tree and one-loop order. 
In the course of study, we have confronted problems of large scale computations 
that prevent the application of the system. The restriction can be removed by the 
hardware development and software engineering. Among them, we have discussed 
here the parallel computation, the singular integral in numerical computation, the 
symbolic manipulation and the floating number with high precision. Armored with 
those items, we can obtain the theoretical predictions for a larger class of physical 
processes in the current and future accelerator projects. 
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